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Abstract. Individual nodes in evolving real-world networks typically experience 
growth and decay — that is, the popularity and influence of individuals peaks 
and then fades. In this paper, we study this phenomenon via an intrinsic nodal 
fitness function and an intuitive aging mechanism. Each node of the network 
is endowed with a fitness which represents its activity. All the nodes have two 
discrete stages: active and inactive. The evolution of the network combines the 
addition of new active nodes randomly connected to existing active ones and the 
deactivation of old active nodes with possibility inversely proportional to their fitnesses. 
We obtain a structured exponential network when the fitness distribution of the 
individuals is homogeneous and a structured scale-free network with heterogeneous 
fitness distributions. Furthermore, we recover two universal scaling laws of the 
clustering coefficient for both cases, C(fc) ^ and C ~ n~^, where k and n refer to 
the node degree and the number of active individuals, respectively. These results offer 
a new simple description of the growth and aging of networks where intrinsic features 
of individual nodes drive their popularity, and hence degree. 
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1. Introduction 

Over the past decade, there has been an explosion of interest in complex networks for 
describing structures and dynamics of complex systems. Despite differences in their 
nature, many networks may be characterized by similar topological properties. For 
instance, real networks display highly clustering than expected from classic random 
graphs [1]. Also, it has been widely observed that node-degree distributions of many 
large networks are heavy tailed [2], e.g., exponential and power-law. To understand how 
these phenomena arise, research devoted to evolving networks has rapidly flourished 
[31 m El El [71 E]. The basic premise is that the network will continue to grow at a 
constant rate and new nodes attach to old ones with some possibility. When the newly 
added nodes connect with equal probability to nodes already present in the network, 
the degree distribution of the nodes of the resulting network is exponential. Whereas for 
newcomers connecting to old ones with linear preference of the node degree, Barabasi 
and Albert (BA) observed a power-law distribution of connectivity [3]. 

In the real world, agents, represented by nodes, always age after growth. For 
instance, in scientific citation networks there is a half-life effect: old papers are rarely 
cited since they are no longer sufficiently topical (or they are more often referenced 
through secondary literature). On the World Wide Web popular web sites (for example, 
search engines) will often loose favor to newer alternatives. To study the effect of this 
phenomenon on network evolution, the BA model has been modified by incorporating 
time dependence in the network P[Tni[IIl[l2lll3l[Tl[T5l[l6l[I71[l8lll9]. Dorogovtsev 
and Mendes studied the case when the connection probability of the new node with an 
old one is not only proportional to the degree k but also to a power of its present age t~" 
(where r is the age of a node) [9j. They found that the network shows scale-free (SF) 
behavior only in the region a < 1. For a > 1, the distribution P{k) is exponential. Yet, 
the gradual aging model show lower clustering than realistic networks. On the contrary, 
Klemm and Eguiluz considered evolving networks based on the memory of nodes and 
proposed a degree- dependent deactivation network model [15] which is highly clustered 
and retains the power-law distribution of the degree (but no consideration of exponential 
degree distributions). 

The aim of this paper is to propose a simpler and more fundamental mechanism to 
build networks, including SF networks, while retaining the positive features of aforesaid 
models. The mechanism we propose simulates realistic networks, and can be understood 
analytically. It has been shown that in some networks, the popularity (or activity) 
of a node is essentially determined by the so-called "fitness" [20l EH [22l |23] which 
is intrinsically related to the role played by each node, such as the innovation of a 
scientific paper or the content of a webpage. This allows us to represent the activities 
of individuals by intrinsic fitnesses and suggest a fitness-driven deactivation approach 
to build structured networks. Compared with topological information, the intrinsic 
fitness provides a more natural and appropriate deactivation criterion for aging of nodes. 
We show that depending on the node-fitness distributions two topologically different 
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networks can emerge, the connectivity distribution following either an exponential or a 
generalized power law. In both cases, the networks are highly clustered. Irrespective 
of the fitness distribution, we observe two scaling laws of the clustering coefficient, 
C{k) ~ and C ~ n~^, where k is the node degree and n corresponds to the number 
of active individuals in the network. Hence, this mechanism offers an explanation for 
the origin and ubiquity of such clustering in real networks. 

2. Model 

Rather than connectivity-dependent deactivation dynamics of the nodes developed in 
[T5| \TQ\ \T7\ [T8| [19], the present deactivation model is based on the individual fitnesses. 
For each node i a fitness Xi > 0, the random number drawn from a given probability 
distribution function p(x), is assigned to measure its popularity or activity. The 
deactivation mechanism is characterized by the transition of a node from the active 
to the inactive state interpreted as a collective forgetting of it. 

The network starts from an initial seed of n nodes, totally connected by undirect 
edges, which are all active. Then at each time step the dynamics runs as follows. 

(i) Add a new node i, which connects to m (m < n) nodes randomly chosen from 
the n active ones. By fc'" we denote the in-degree of a node, i.e., the number of edges 
pointing to it. The in-degree of the newcomer is kf^ = at first. Each selected active 
node j receives exactly one incoming edge, thereby A;™ — )■ A;™ + 1. Since the out-degree 
of each node is m always, the total degree of a node is k = /c™ + m. 

(ii) Activate the new node and deactivate one (denoted by j) of the n old active 
nodes with probability 

TT{xj) = axj^, (1) 

where a = iJ2jeAXj^)~^ is the normalization factor. The summation runs over the set 
A of the n old active nodes. 

During evolution, a node might receive edges while it is active, and once inactive it 
will not receive edges any more. Note that the fitter the individual is, the more difficult 
for it to be deactivated. For the case of the citation network, Eq. ([1]) means that the 
famous paper with great innovation is less possibility to be forgotten. 

3. Degree distribution 

Denoting by a{k, x, t) the probability of active nodes with degree k and fitness x at time 
t, we can write out the master equation for network evolution 

m m 
a{k + l,x,t+l) = a{k + l,x,t)[l-7i{x)]{l-—)+a{k,x,t)[l-'K{x)] — .{2) 

At each time step, an active node is deactivated and a newcomer joins the set A to keep 
the number n unchanged. According to Eq. ([1]), the normalization factor a varies with 
time because of the change of the nodes in the active set. However, the fitness of each 
node is a random number taken from a given probability distribution p{x), therefore the 
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normalization a fluctuates very slightly and can be treated as a constant. Substituting 
Eq. ^ into Eq. ^ yields 



a{k + l,x,t + l) ^ a{k + )(1 



X 



—) + a{k,x,t){l--) — .{3) 
n X n 



Imposing the stationarity condition a{k,x,t) = a{k,x), we obtain the equation 

H k—m 



a{k + 1, x) 



mix — a] 



mix — o" 



-a{k, x) 



mix — a] 



mix — a) + na 



a{m,x) (4) 



for the probability of active nodes with degree k + 1 and fitness x in the stationary state, 
where 



a(m, x) 



nix — a) 



-p{x) 



(5) 



m{x — a) + na 

is the stationary probability of active nodes with degree m and fitness x. Then we 
obtain 

^ k~m 



a{k, x) 



n 
m 



mix — a) 



mix — a] 



na 



p{x). 



(6) 



Denoting by a{k) the possibility of active nodes with degree k in the steady state, we 
have 

k—m 



a{k) 



a{k, x)dx = — 
Jo 



m 



m{x — a) 
mix — a) + na 



p(x)dx. 



(7) 



In case the total number of nodes in the network is larger than the number n of active 
nodes, the degree distribution P{k) can be approximated by considering inactive nodes 
only. Thus, P{k) can be calculated as the rate of the change of a(/c), 

daijx) 



P{k) 



dk 
m 



mix — a] 



k—m 



In 



mi^x — a) 
mix — a) + na 



p{x)dx. 



m{x — cr) + na 

Even when the form of p{x) is given, it is still difficult to solve the integral on the 
right-hand side of the equation. Instead, we need a more subtle technique. We assume 
that 

m{x — a) 



F{x) 



mi^x — a) + na^ 



n 

G{x) = -[F{x)]-"'ln[F{x)]~'p{x) 



m 



(9) 

(10) 



Without lack of generality, we also normalize the fitnesses. Now Eq. ([8]) can be rewritten 



as 



P(k) = / [F{x)]''G{x)dx. 



(11) 

As will be seen below, for the proper choice of F and G, one can construct networks 
with exponential or power-law degree distributions, and then determine the forms of 
the corresponding fitness distributions. 
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(i) Exponential degree distribution. We set F{x) = l//i and G{x) = u, where fi{> 1) 
and u are positive constants. Consequently, the integral of Eq. ffTTl) is 

p{k) = z//i-^ (12) 

following an exponential. According to the definition of F{x), we have 
mix — a) 1 



m{x — a) + na fj,^ 
the solution of which reads 



(13) 



na 

X = a -\ 7 r = constant. (14) 

m(/i — 1) 



Then the normalization factor becomes 

-1 



a = <n 



na 

a 



m{j^ — 1) 



-1 



(15) 



yielding 



m{n — 1 



According to the definition of G{x), we have 
which suggests 



n 

u = —fi"'\nfip(x), (17) 
m 



mu 

P(^) = '-i = constant. 

n/i'" In /i 

According to the normalization, p{x) should satisfy Jq p{x)dx = 1. Combining Eqs. ([1 
and (ITS]) , we finally obtain the distribution of p{x), 

^( a(ma — m + n)\ , , 

^w-n^" 4-1) ); 

That is to say, the fitness of each node is identical. Considering the constant 
approximation of a, the above result can be generalized to homogeneous distributions 
of fitnesses. Here, homogeneity implies that the node fitnesses have small fluctuations 
around the mean (x), hence the small variance. We conclude that for any given 
p{x) distributing homogeneously, the fitness-driven deactivation generates structured 
exponential networks. 

(ii) Power-law degree distribution. We set F{x) = e~'^^ and G{x) = ipx'^, where 0, 
(/? and 7 are positive constants. Accordingly, Eq. (fTTl) can be rewritten as 



Pik) = 7-^l e-z^dz. (20) 



(t>k 



In the limit k — )■ oo, one has J^^ e~^z'^dz = r(7 + 1). For any finite k, the integral 
is convergent and smaller than r(7 + l). Thus the connectivity distribution 
has a power- law form 

P{k) ~ k-^-\ (21) 
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Figure 1. (Color online) Degree distributions of nodes of generated networks for 
different fitness distribution functions: uniform (a), Gaussian (b), exponential (c), 
and power-law (d). All the fitnesses have been normalized. For uniform and 
Gaussian distributions with the same mean (0.5) and small variances (0.08(4) and 
0.01(6)), the node fitnesses fluctuate slightly around the mean and can be regarded 
as the homogeneity; whereas for exponential and power-law distributions, there are 
large fluctuations of the node fltnesses due to the right-skewed feature, resulting in the 
heterogeneity. The experiment networks start from the initial n = 10 nodes and end 
with the total population N — 10^. 



Combining Eqs. ([9]) and ( fTOl) . we obtain 

(f)n 

This heavy tailed distribution imphes large fluctuations of the node fitnesses, hence 
heterogeneity. We conclude that the fitness-driven deactivation can build structured SF 
networks if the node fitnesses are distributed heterogeneously. 

In Fig.[T]we give the simulation results of degree distributions P{k) for four kinds of 
distribution functions p{x) of node fitnesses. In case that the distributions of fitnesses 
are homogeneous, e.g., uniform (Fig. [H^a)) and Gaussian (Fig. [11(b)), the linear-log 
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Figure 2. (Color online) The dependence of the slopes of the lines in Fig. [T] on m. In 

the inset, the data reported on the log-log representation show a power law m~^, and 
the best linear fit gives ^ = 0.98(1). 



plots imply an exponential degree distribution. Conversely, when the fitnesses distribute 
heterogeneously, the log-log plots of Figs. [H^c) and[T](d), corresponding to exponential 
and power-law cases respectively, predict a generalized power law. We obtain the slopes 
of the lines in Fig. [1] by least squares fitting and plot them as a function of m in Fig. [2j 
For homogeneous fitnesses, we notice a scaling relation between the slope and m, whereas 
for heterogenous fitnesses, the slope shows independence on m. 



4. Clustering coefficient 

In a network, if a node i has fcj edges, and among its ki nearest neighbors there are Cj 
edges, then the clustering coefficient of i is defined by 

In the deactivation model, new edges are created between the selected active nodes and 
the added one. Let us first consider the case of n = m. At each time step, the degree 
ki of the active node i increases by 1 and Cj increase by m — 1 until it is deactivated. 
Therefore, the evolutionary dynamics of ki and Cj are given by 

ki ={m + t), (24) 

(25) 
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Integrating Eq. fl25p with the boundary condition ej(0) = m(m — l)/2 and substituting 
the solution into Eq. (|23|) . we recover the clustering coefficient c{k) restricted to the 
nodes of degree k [161 E] 

2(m — 1) m(m — 1) , , 

For n > m, the clustering coefficient C{k) is just the generalization of Eq. (126|) . 
_ r2(m-l) m(m- 1) ] 

n [ " k{k-l)\' ^^^^ 

which indicates that the local clustering scales as C{k) ~ k^^ for large k. The clustering 
coefficient C of the whole network is the mean of C{k) with respect to the degree 
distribution P{k), 



C= / C{k)P{k)dk. (28) 

J m 

Substituting Eqs. (fT2l) and (1211) into the integral respectively, we obtain 



{ 2i/(m^-m) roo fi , u{m'^-m^) roo ^ 

n Jn^ k n ■J'm k{k—l) ' (OQ^ 

2{'m?-m) roo dk , jm^-m'^) roo dk ^ ' 

n Jm fc7+7 "T „ Jm fc7+2(fc_i) i 

proportional to n^^. 

In Fig. [3] we show the simulation results of C{k) as a function of k. All the plots 
follow a power law C(A;) ~ k~^, which coincides with the expression in Eq. f l27|) . For 
the clustering coefficient of the whole network, as shown in Fig. HI the linearity of all 
the plots also implies a power-law relation C ~ n~^. It is worth noting that both the 
scaling laws are independent of fitness distribution functions. One obtains the same 
result for uniform, Gaussian, exponential and power-law distributions of ffinesses. 



5. Conclusion 



In this paper, we have presented an alternatively simple and intuitive model for a large 
and important class of networks widely observed in the real world. We defined the 
intrinsic fitness as the way of quantifying the popularity of individuals, i.e., the ffiter 
a node is, the higher possibility it is active. The growth dynamics of the network 
is governed by the naive fitness-driven deactivation mechanism. The deactivation 
probability of a node is proportional to the inverse of its ffiness, which characterizes the 
individual capability of obtaining further links. We studied the connectivity distribution 
and the clustering coefficient that can fundamentally shape a network. On one hand, 
we found the great influence of the node ffinesses on the connectivity distribution. The 
homogeneous ffinesses generate exponential networks, while the heterogeneous ffinesses 
result in SF ones. On the other hand, we recovered two universal scaling laws of the 
clustering coefficient regardless of the ffiness distributions. These results are consistent 
with what has been empirically observed in many real-world networks in (for example) 
[211 [251 ISSl [27] , and so the present model provides a new way to understand complex 
networks with age. 
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Figure 3. (Color online) Log-log plots of C{k) as a function of k for different fitness 
distribution functions: uniform (a), Gaussian (b), exponential (c), and power-law 
(d). The solid lines correspond to a power law 2k~^. 
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